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Abstract. We prove involutivity of Einstein, Einstein-Maxwell and other field equations by 
calculating the Spencer cohomology of these systems. Relation with Cartan method is traced 
in details. Basic implications through Cartan-Kahler theory are derived. 
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^ ' Introduction 
> I 

O , The field equations, derived independently by A.Einstein [E] and D.Hilbert [H], form the basis 

' of the general relativity. Soon after these equations had appeared in 1915 (see |LMP1 IT] for the 

iy»^ , history of invention) some very important particular solutions were found, yet not much has 

been known on their solution space before E.Cartan's contribution. 

In correspondence with A.Einstein [CE| he established involutivity of various field equations. 
Cartan approached this through his theory of exterior differential systems |Ci| , namely by calcu- 
lating Cartan characters and verifying that they pass the Cartan test. This is quite an involved 
work. Detailed proof of this is contained in [Ga| (see also [DTj : quite a different proof of invo- 
(— I , lutivity for Einstein, Yang-Mills and other equations is contained in |DV| ) . 

In this paper we prove involutivity of the field equations using the formal theory of differential 
equations [S] : we calculate all Spencer cohomology of the system and check their vanishing in the 
prescribed range (together with vanishing of the structure tensor). By Serre's contribution to 
|GS| this is equivalent to Cartan test. Since chasing diagrams is considered nowadays standard, 
' this turns out to be a reasonable path. 

^ , Let us notice that up to now the Spencer (5-cohomology has never been evaluated for the 

field equations. We do our calculations for both Einstein and Einstein-Maxwell system, and 
\ also encompass such systems as pure radiation and dust. We start with vacuum or source- free 

equations and then impose additional fields, proving involutivity of various field equations. 
We relate our calculations to those of Cartan which is not obvious, since the two theories 

■ - Cartan and Spencer - though accepted being equivalent, are not in direct correspondence. 
I Finally we derive some simple but important implications using the Cartan-Kahler theorem. 

■ Let us mention that in his papers IC2I Cartan mostly considers the so-called unified field theory 
based on distant parallelisnu, which corresponds to Einstein system with non-zero torsion, so 
that the number of differential equations in the system is 22, not 10 or 18 (or 14) as in the usual 

^ . , Einstein and Einstein-Maxwell equations. Involutivity of these two latter do not follow from 

' involutivity of the former upon a specification. In addition, Cartan arguments by exhibiting 

relations between the equations but not proving they are all. The formal theory approach, 
adapted here, provides both rigorous and economic way to prove involutivity. 

1. Background: jets. Spencer cohomology and all this 

We will consider here only the theory of systems of PDEs of the same order k |Go[ [S] . The 
general theory, developed in |KLi[ [KL3 I, shall be useful for other purposes. 

Thus let E — Ek C J'^tt be a submanifold in the space of fc-jets of sections of a bundle tt : E 
M, subject to certain regularity assumptions, which include the claim that Hk,k-i ■ £k ~* J^~^tt 
is submersion. We let = J'tt for I < k and £1 — sj! ''^ for I > k, where the latter space is the 
prolongation defined as 

^k~^^ — {Mx G J^''^ ■ jet-prolongation j''{s) is tangent to £k at [s]^ with order (/ — fc)}, 

where [s]^ is the k-jet of the (local) section s G r(7r). Equation £ is called formally integrable if 
all the projections iri^i-i : £1 — > £i-i are submersions. 



Key words and phrases. Einstein equations, Einstein-Maxwell equations, involutivity, Cartan numbers, sym- 
bols, Spencer cohomology MSG: 83C05, 83C22, 58H10, 58A15. 

^Specific references are vol.11 p. 1199-1229 and vol.III-1 p.549-611. 
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Let US denote by N the tangent space to the fiber of tt and by T the tangent space to M. 
Then the symbol spaces gt C S^T* N are the kernels of dTTt,t~i ■ TEt TSt-i- We obviously 
have gt = S^T* (g) A'' for t < k and the space gk is determined by the equation, however the 
higher index spaces are difficult to calculate without knowledge of formal integrability. 

Instead one considers formal prolongations defined as gt — {gk ® S*'~''T*) n (S'*T* ® N) for 
t > k. These symbols are united into the Spencer (5-complex 

O^gt^ gt-i ®T* gt-2 ® A^T* ■ ■ ■ gt_, ® K'T* ^ . . . (1) 

with morphisms 5 being the symbols of the de Rham operator. The cohomology at the term 
gt-i (8> A*T* is denoted by H*~^''^{8) and is called Spencer (5-cohomology of £. 

Formal theory of PDEs describes obstructions to formal integrability as elements Wt S 
H^~^'^{£), called curvature, torsion, structure functions or Weyl tensors. Their vanishing is 
equivalent to formal integrability (and in certain cases to local integrability). 

Symbolic system g = ®gt is called involutive if H^'\g) = for alH ^ fc — 1 and i + j > 0. 
This is equivalent to fulfillment of Cartan test for the corresponding EDS (which in turn means 
a PDE system of the 1st order). 

Equation £ is called involutive if its symbolic system is involutive and in addition the only 
obstruction Wk vanishes. Thus involutive systems are formally integrable. 

Advantage of involutive systems is that compatibility conditions should be calculated only 
at one order, while in general they exist in different places and one shall carry the whole 
prolongation-projection method through [Si IKLVl IKL3I . Fortunately many equations of mathe- 
matical physics are involutive (this can be easily checked for all determined and underdetermined 
equationtp) and we are going to prove this for relativity equations. 

2. Einstein equations 

We run the setup very briefly, referring to plentiful books on differential geometry and rela- 
tivity for details, e.g. PIHISt]. 

Let M be a (4-dimensional in physical applications) manifold, g pseudo-Riemannian tensor 
(for relativity: of Lorentz signature (1,3)) with Ricci tensor Ric and scalar curvature i?, A a 
cosmological constant and T the energy-momenta tensor. The Einstein equations [El |H] are: 

Ric-\Rg + Kg = T. (2) 

We will assume at first that T is a given tensor (in some models like electromagnetic it is 
traceless, which implies that the scalar curvature R is constant, but in general it is not so), and 
only the metric g is unknown (further on we'll treat the case, when T depends on unknowns 
fields entering the equations). We can re- write ([2]) as 

G(Ric) = Ta, 

where G(h) — h — ^ Tig{h) g, h € S^T*, is the gravitational operatoiH and T\ = T — Ag. 
Bianchi identity reads 

t)(Ric) = 0, where d = div^ o G (3) 

and diVg is the divergence operator on symmetric tensors, so that Dh{S_) = Trg(V./i(-, ^)) — 
iV^(Trg h) for £, € T. This implies the following conservation law: 

diVg(rA)=diVg(r)=0. (4) 

This is the first order PDE w.r.t. g and so system ([2]) is not involutive unless T = (indeed, we 
get a non-trivial equation of lower order than ([2]), see Appendix for more details on this). 

Thus in what follows in this section we'll concentrate on the vacuum cas^: T = 0. 

Tracing ^ by g yields 4A = i? = const, so that the Einstein equation £ is equivalent to 

Ric = A.g. (5) 



An equation is (under)determined if the symbol of its defining operator is an isomorphism (resp. surjective). 

^Herc and below T = TM is the tangent bundle to M and T* is the cotangent bundle. S^T* C S'^T* is the 
bundle of non-degenerate quadrics (not necessary positive definite). 

^DeTurck's idea [B] is to use covariance of the left hand side G[g] of ((2]l and change the equation to G[g] = </3*T, 
where (/? : M — » M is a diffeomorphism, so that T is given while [g, (p) unknown. This system (coupled with 
compatibility S^p gT = 0) is already involutive for any non-degenerate T - see Section [5] DeTruck |DTI proves 
solvability of ^ differently - by establishing involutivity of its prolongation-projection l2ll+ll4t. 
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To understand this equation we need to study solvability of the Ricci operator 

Ric : C°°{SlT*) C°°{S^T*), 

which gives rise to the sequence of operators 4>B.ic ■ J^'^'^ {S\T*) J^{S'^T*) with symbols 

c^Ric : S''+^T* (g) S^T* S^T* ® S^T* . The symbol of the Ricci operator was described in 
|DT[|Ga[[B] and it equals (we choose p,q, . . . £ T*; it suffices to define the map on decomposable 
quadrics) 

CTRic(/ <E> q^) = {p,q)gpq - ^dplg?^ + kl^), 
which is clearly an epimorphic map for any signature of g. The prolongation 

41ikT2P'^^ ® = {p,q)gp'r^pq ~ U\P\>' ® 9' + klg/ 

is not epimorphic for fe > as follows from the Bianchi identity. 
Symbols of the Einstein equations (O are precisely 

gk+2 - Ker(4';^), fc > 0, 

and we let 50 = S^T* , g-^=T* ® S^T*. 

We calculate the Spencer cohomology of £ ©by constructing resolutions to the symbols of 
the Ricci operator. The first Spencer complex is exact. The second Spencer complex is included 
into the commutative diagram, implying H^-^{£) — S^T* and ) = for i 7^ 1: 





4 I 4- 

> 92 >■ gi (g) T* y g„ « A^T* v 

S^T* ® S^T* T* ® S^T* (g T* ^ S^T* A^T* 

^ S^T* > 



In what follows we shorten S'^T* to S'^, and use similar notation A'^ — A'^T* for readability of 
the diagrams. The third Spencer complex is included into the commutative diagram, implying 
H^^^{£) = T* and H^-'''{£) = for i 7^ 2: 





1 i i i 

>■ 93 >■ 52 «> T* >■ A2 > go ® ^ 

-> T* 5*2 52 T* > 

(Ti, I ^ / I 

y T* V 



The formula for the symbol of operator © d : C°°{S'^T*) C°°{T*), 
a:,{p(^q-r) = ^{{p,q)gr + {p,r)g q - (q,r)gp) 



easily implies exactness of the first column. The next commutative diagram is exact: 
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\^ ^ ^ 

^ .94 ^ 53 (g) T* ^ 52 ® ^ 51 (g) ^ 50 (g) A'^ ^ 

^ 5*4 5*2 y 5*3 (8) 52 (g) T* ^ 52 ® 52 ® A2 ^ T* (g) S'2 (g) A^ — v 52 ^ ^4 ^ q 

'^Ric'l' i I ^ 

^ 5*2 g) 52 — ^ T* ® 5*2 g) T* — >■ 5*2 g) A2 > 

T* ®T* > T* g) T* > 



and this extends to the commutative diagram for any fc > 0, with exact rows and columns (we 
draw the diagram for the case n — A, but this has an obvious extension): 



^5fe+4 ^5fc+3(gr* ^5fe+2®A2 ^5fc+ig)A3 >g^(g,K^ > 

^ S"'-+4 5*2 ^ S"^'+3 g) 5*2 g) T* ^ S"'-+2 5*2 g) A2 ^ S'-^+i g) 52 g) A^ ^ 5"= g) 52 g) A" ^ 

o-^j^ 'I I ^ \ 1 

^S'''-+2 5'2^5''=+ig)5'2g)T*^5''=g)5'2g)A2^5'''-ig)5'2g)A3^5'=-2®S'2 0A'^^ 
^S''=+i ®T*^S''= g)r* g)r*^S''=-i ® A2^S''=-2 0T* ® A3^S''=-3 g)r* g) A*^ 



Since H^'^{E) = for fc > 2 and aU I, the symboHc system 5 = ©5fc of f is involutive. To 
prove involutivity of PDE system £ it is thus enough to check that the actual compatibility 
conditions belonging to H^'^{£) vanish. But these elements, as follows from calculation of their 
symbols, coincide with Bianchi relations and thus their equality to zero holds identically. 

Notice that we don't use specific features of Lorentz geometry, and our arguments work 
generally (though the diagrams become larger). We have proved: 

Theorem 1. The only nonzero Spencer 5-cohomology 0/ Einstein equation (0)t=o 0*^?/ di- 
mension and signature are 

H"'%£) ~ S^T*, H^^\£) ~ S^T*, H^^^S) ~ T*. 

The Einstein equation £ is involutive. 

3. Relation with Cartan approach 

In Cartan theory involutivity is checked via Cartan characters, which are defined as follows. 
Consider a symbolic system 5 of order k: 

5,; = S'T* ®N, i<k; 5fc C S'^T* g) iV; 5, = 5^'-'=^ C S'T* g) TV, i > k. 

Let — V* C V*_i C ■ ■ ■ C Vq = T* he a generic complete flag (difference of dimensions is 1; 
n — dimT = diniAf). By definition 

= dim(5fc n S''V*_^ (g> N) ~ dim(5fe n S'^V* (g, N), l<i<n 

(the sequence monotoni cally de creases) , so that dim5fc = si + • • • + s„. 

Cartan test |C2[ IMai iBCG'^j claims that symbolic system 5 is involutive iff 

dim5fe+i = si + 2s2 H h ns„. 
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In this case we also have dimgfe+i — s'l + ■ ■ ■ + s'„, where — Si + ■ • ■ + s„ are the Cartan 
characters for the prolongation, i.e. 

s', ^ dini(.gfc+i n S''+^V*_i ® iV) - dim(.gfe+i n S'^+'^V* ® iV), 1 < i < n. 

Thus we can calculate the dimensions of symbol spaces via Cartan characters: 

a».,.„<'-ro.-<. .... !:('-::;->„... 

where m = dimA'^, which we also denote by sq. 

Let us relate Cartan characters and Spencer J-cohomology for an involutive system of pure or- 
der k. The only nontrivial dimensions of the latter are ho = dimiJ''^" = m and hi = dimiJ'''^^^*, 
1 < i < n. To one side the relation is given by 

Proposition 1. The numbers {ho, . . . , ft,„) are expressed through (sq, ■ ■ • , Sn) cis: ho = sq and 

'.,=(-iyE..|:(-i)'(:)('+5:r>(-i)'-D-i)<:)C'-':!r'0 

for 1 < i < n. 

Proof. Due to involutivity the Euler characteristic of Spencer complex ([T]) equals {—ly^'^'^-^ ht-k+i 
for t > k, zero for < t < A: and ho for t = 0. Calculating it directly as X]r=o(~l)* dim(7t_i(") 
and using ([6]) we get the result. □ 

The relations above are invertible, but we obtain the inverse formula from another idea. 

Proposition 2. The numbers (sq, . . . , Sn) are expressed through (ho, . . . , /i„) in triangular way: 
So = ^0 o.iT-d 



si = 



^ i=n-l + l ^ ' 



Proof. For I ^ 1 the expression H^iV) — 'Ylii<i dimg,; is a polynomial, called Hilbert polynomial 
of g and so dim^^ = Hgi^z) — Pls(z — 1) is a polynomial too (for large integers z = I and we 
extend it to the space of all z e C). 

The Hilbert polynomial can be computed through the standard resolution of the symbolic 
module g* ^Grl IKL2I and we get: 



dim = Y^i-'^y dim H''\g) 



z + n — i— j — 1 
n — 1 



n-l n-1 )+H n-1 l-"' 



On the other hand from ^ we have the following expression: 

dim,.=E^.(^+;:^-^). (8) 

Comparing ([7|) to ((S]) we obtain the result: At first substitute z = k — 1 and gelQ 

fn + k-2\ 
si = fto ] - h„ 



then calculate difference derivative by z, substitute z = k — 2 and get the formula 

/n + fc-3\ , , fn-l 

S2 = ftO „ - hn-l +hn[ 

\ n — 2 J \n — 2y 

and so on. □ 

Remark 1. To see that formulae of proposition\^ invert these of proposition{l\is not completely 
trivial: one must use certain combinatorial identities. 



^One shall be careful: in this substitution dim gz is understood as analytic continuation because the actual 
value of dim(;j._j could be different; on the other hand studying the large values / one gets the same result. 
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Now let US apply the result to Einstein vacuum equations (we restrict to the physical dimension 
n = 4, but due to previous formulae the general case is easily restored). As we calculated in the 
previous section 

ho = 10, hi = 10, /i2 =4, /i3 /i4 = 0. 
Thus proposition [5] implies that the Cartan characters are 

si = 40, S2 = 30, S3 = 16, S4 = 4. 

This calculation was independently verified in Maple with (Differentia I Geometry): [X] , 

In particular, the Cartan genre is 4 and the Cartan integer is 4, i.e. the general (analytic) 

solution of the Einstein vacuum equations depends on 4 functions of 4 arguments. This is indeed 

so due to covariance: the group Diffioc(M) acts on £ as symmetries. 
We can calculate the Hilbert polynomial of the Einstein equation 

6 6 

The first dimensions of the symbol spaces are: 

dim go 10, dimgi = 40, dim 32 = 90, dim gs — 164, dim 34 = 266, dim 175 = 400, . . . 
(this in particular shows that direct calculation can be costly). The Cartan test works as follows: 

si + 2s2 + 3s3 + 4s4 — 164 — dim 53. 

Historical remark. Explicit formulae for Hilbert function and polynomial of PDE systems 
were calculated already by Janet |jj (and maybe to some extent were known to Riquier). 

Nowadays his result is generalized to arbitrary involutive bases and this provides a link be- 
tween Cartan characters and Spencer cohomology. The first step - to read off the Hilbert function 
from an involutive basis (or Pommaret basis or Cartan characters) - can be found in a number 
of papers [Si [Sg [GBl El] . 

The second step comes from commutative algebra, where it is well-known that the Koszul 
homology and the Hilbert polynomial are related (beware that involutive systems are more 
important in PDE context and are infrequently discussed in algebraic situations). Since Spencer 
cohomology are M-dual to Koszul homology, this yields a principal link, mentioned above. 

This link however is not clearly marked in the literature, and the above explicit formulae are 
new. In addition we deduce relations in both direction without implicit inversion (involving lots 
of combinatorics). 

4. Einstein-Maxwell equations 

These equation extends ^ in the sense that the energy-momenta tensor is prescribed as 
electromagnetic tensor. Denote by J the current density. Einstein-Maxwell equations have the 
following forrr0: 

Ric -ii?.9 = {F^)o, dF = 0, 6gF = J. (9) 
Here the tensor F in the first equation is viewed as a (l,l)-tensor (an operator field) via the 
metric, and (-F'^)o — F^ — -j Tr(F^) is the traceless part of its square, while F in the latter 
equations is a 2-form, and 6g — ±* d* : it^M fl^M is the Hodge codifferentiafl. 

In order not to deal with involutivity of systems of PDEs of different orders (the theory 
developed in |KLi| ), we can re- write the system as a pure 2nd order system by introducing the 
potential A G n^M, F = dA: 

Ric -ii?.9 = {dA o dA)o, Sg^dA) = J. (10) 

Both systems ([9|) and (fTO|) have the following compatibility condition of order 1 in 5: 5gJ — 0. 
Thus they are not involutive unless J = 0. This will be assumed at the end of this section. 

But let us study at first the pure Maxwell equation (with known g), written as a 2nd order 
system with operator II = Sg o d: 

H(A) = J. (11) 

The symbol of this operator equals 

an-. S^T* (»T* ^T*, Q ® p ^ p^Q ~ TTg{Q)p, 

^We set the cosmological constant A = 0, which docs not restrict mathematics (can be incorporated back 
without destroying any conclusion), but agrees with physical observations. 
^Not to be confused with Spencer (5-differential. 
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where in the first term to the right the duahzation T* ~ T is used. Thus the symbol is epimor- 
phic, while its prolongations are not, since they have left divisor of zero: 

f^^r^^ ° ^ = ^^"^ '^s, : T* ® T* ^ M, q (g, p ^ {p, q)g. 

The symbol of 5g is however epimorphic together with all its prolongations and so we get the 
sequence of commutative diagrams with all rows and columns exact except for the top (Spencer 
5-complex) and the bottom rows: 



\ \f ]r 

>■ 92 > .91 ® T* > g„ (g) A^T* ^ 

>■ (g) T* ^ ^ > T* g) T* ® T* y T* A2 ^ 

'^n I ^ - " i I 

^ T* ^ 

I 


This imphes iJ^'^ ~ T*. The next complex 


\- \- ^ \- 
> .93 y .92 T* y gi ® y ® A^ ^ 

^ S'3®r* — y S"^ ®T* ®T* ^T* ®T* ® y T*®A3 

^ T* T* ;r^^ T* ® T* ^ 

'^^s I ^ / .|, 
y M 

\ 


yields iJ^'^ ~ M. Further complexes are already exact. Here's the next one: 



1^ 

y .94 y 33 ® T* y g^ ® A^ y .91 (g) A^ y go A'^ ^ 

^ S^^T* — y 5*3 (g) T* (g) T* ^ 5*2 ® T* A2 ^ T* (g) T* (g) A3 — y T* ® A^ ^ 

^ 52 T* — ^ T* g) T* g) T* ^ T* ® A2 ^ 

y T* y T* y 



and one can easily prolong. Since dimiJ*'^ = 1, there is only one compatibility condition and 
from its symbol we identify it with the condition 5gJ — Q (which comes from the Hodge identity 
5^ =0). Since the metric g is fixed this implies: 

Theorem 2. The only nonzero Spencer 5-cohomology of Maxwell equation f— n in any di- 
mension and signature are 

The Maxwell equation is involutive provided the compatibility 5gJ — holds. 
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Now we can study Einstein-Maxwell equation. Energy- momentum tensor T = (^■^)o satisfies 

divgr = F(.,j«), 

where j" is the g-Mt of J e fi^Af ([K ). Therefore compatibility conditions for equation ^ or 
(|10p are: j") = and 6gJ = 0. They are trivial only for J = 0, and since they have lower 
order (in g and F) this must hold for involutivity. 

Remark 2. It is tempted to change the Einstein-Maxwell system in the case J ^ 0, similar to 
DeTurck tric^ to the following: 

Rfo-ii?g = (F2)o, dF^O, 6gF ^ (j)*J 

for the unknown {g,F,(j)). While the second compatibility condition Sg{(j)*J) = yields an un- 
derdetermined equation of 2nd order on cf), the first condition F[-,(f)*.P) = implies the 0-order 
condition dot F = Q on F, and so the system cannot be involutive. 

Thus we should restrict to the case of no external sources: J = 0. The key observation is 
that this system is weakly uncoupled, meaning that its symbol splits into the sum of symbols 
of Einstein and Maxwell equations. Thus Spencer cohomology becomes the direct sum, and the 
compatibility condition for the system £M ()10p is the union of two respective compatibility 
conditions (Bianchi and Hodge identities). 

Theorem 3. The only nonzero Spencer S- cohomology of source- free Einstein-Maxwell equation 
ilO\) r-n in any dimension and signature are 

H"^°{£M) ~ S^T* ® T*, H^^^iSM) ~ S'^T* © T*, H^''^{EM) ~ T* © R. 

This Einstein- Maxwell equation EM is involutive. 

Couple of remark to properly place this results are of order. 

Remark 3. With Einstein-Maxwell system (0) we get H°'" = .go ^ 5^ © T* ® T* , 

which correspond to Rainich "already unified field theory" [R] ('so that we don't have bosonic or 
fermionic parts, but just tensors). The other cohomology do not sum (since belong to different 
bi-grades) but unite and we get that the only non-vanishing Spencer cohomology of 0) are: 

~ M© A^'T*, H^'^ ~ T*, 

^0,3 ^ ^brji* F[^'^ ^ ^6^^* 

Notice that for n = 4 the latter line disappears. 

Remark 4. Following Rainich [11 Einstein-Maxwell equations are equivalent to the system 

(Ric2)o = 0, i? = 0, 

where Ric is viewed as operator and Lq ~ L — ^Tr{L) is the traceless part of an operator L. 
Though relation between the two systems is non-local (but rather simple integral), involutivity 
holds for them simultaneously (however the Spencer cohomology vary, as well as Cartan charac- 
ters). The latter system, though more compact, is fully non-linear and is more complicated. 

5. Relativity equations of other types 

Theorem 4. Consider the following PDE system £ 

G[Ric[g)]^T, L[g,4>]^Q, (12) 

where T = T[g, 0] is an operator of order 1 in the fields (j) and the gravitation (metric) g, and 
L : C°°{M,Ni) — > C°°{M;N2) is an (under) determined operator of order r in (p when g fixed 
(dim A^i > dim Af2/ it is also a differential operator in g of order < r; we assume r — 1,2). 

Suppose that the Bianchi identity diVg(T) = is a differential corollary of equations L — 0. 
Then the system £ is involutive and its only non-trivial Spencer cohomology groups are 

H°'"i£) = S^T* ® Ni, H*'\£) = S^T* ® N2, H^^^{£)^T*. 



'It contains 10 equations on 10 unknowns, the same as for Einstein equation l(2]|. 
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For r = 2 we have H*'^{£) = H'^'^. For r = 1 we get H*^^{£) = H^'^®H°'^ and the system is 
of mixed orders in the sense of |KLi| . It is possible to treat the operator L of mixed (and high) 
orders ri < ■ ■ ■ < using the same technique, but we omit this for transparency of exposition. 

Let us notice that the theorem assumption on compatibihty is equivalent to diVg(T) e [[L]], 
where the latter denotes the differential ideal of L. This happens due to (under)determinacy of 
L and implies that the symbol crjivg t is divisible by ■ 

Proof. Let us restrict for shortness of the diagram^ to the case r = 2 . Similar to Section we 
derive the claim from the bi-complex 



1 i i ^ 

>■ 9k+3 >■ gk+2 «) T* y g^+i y ® A^ i 

\- \- \- \- 

-y S''=+3 (g)Vi^ 3''+^ (g,Vi(g>T*^ S^+^ ® Fi ® A^ ^ 5^= ® Vi A^ — v • ■ • 

"'gIbl i i i i 

S-'^+i ® ^2 ^ 5''= (8) V2 (8) T* ^ S''-^ ® F2 (8 A2 ^ S''-^ (g>V2 (S ^ ■ ■ ■ 

'^dtveo \- \- I 

— > 5*^ 8) T* ^ S''-^ ® T* T* ^ S''-'^ ® T* ® A2 ^ S''-^ ® T* A^ -> • ■ ■ 

i 4- ~l' \- 

00 

where Vi — (B Ni and V2 = 5^ © A^2- The columns are exact and all the rows except the top 
one are exact as well (for all k but —1 and 0). This claim has to be checked only for the fist 
column at the place 5''^+^ (8 V2. 

The prolonged symbol of the equation has triangular form (we omit the sign of prolongation) 

IG 

q d) 



and the next operator is 



(Tdiv 





(S-'^'+i 5*2) © (5*^+1 ® N2) {S'' ® T*) © 0. 



Suppose (Tdiveo(^i: ^2) — 0, i.e. div(6i) = 0. We wish to solve 17(30^(01,02) — (61, ^2) that is 
equivalent to adai) — bi and CT^(a2) = 62 — f£(ai). The first equation is solvable by the results 
of Section [2l The second is solvable by the assumptions of the Theorem. 

The claim about Spencer cohomology follows and we conclude symbolic involutivity. The 
compatibility condition in H^''^{£) is precisely the same as before - Bianchi condition ~ and it 
holds identically due to the assumptions. □ 

Let us now specify certain important sub-cases of the above theorem. In some occasions we 
will need to assume the Lorentzian signature (more generally: non-definite 17) and any dimension 
n > 2, while for the rest the metric tensor can be arbitrary. We refer to [Kj [St] for physical 
motivations. 

1. Pure radiation (null dust). In this case T = ek^ , where fc is a null vector field 
= 0, fc'' = (fc, ■)g its (7-dual and e > a scalar function (it can be absorbed into k, but we 
refrain from doing it). 

The compatibility divg T = is equivalent to the condition 

Vfcfc — Xk, where A — —Ck{t) — edivg(fc) 

{Cy is the Lie derivative along v). Thus fc is a null pre-geodesic and properly choosing e (re- 
parametrization) we get A = and the above under-determined equation becomes determined: 
the first order operator L[g^k,f\ splits into ii © L2 with 

Li = Vfcfc = 0, L2 = divg(fc) + £fe(loge) = 0. 



^For other r one shall change 5*^+^ ® Vi to S^^'^ ® S'^ ® r ^^^^ have similar changes elsewhere. 
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Denoting by £ ^ aL : T* (g) {T i 
on a covector p: 

ip{n, e) 



r©R the symbol of this operator at (fc, e) we get its value 
8 (K,e)) = i{p,k)K, {p,k) + {p,k)^). 



\u\l 



Thus det[^p] = (p, fc)"/e is non-zero for p ^ Ann(fc) and the operator L is indeed determined. 
As an effect, the equations of pure radiation are involutive. 

2. Perfect fluid. Here T — {e + P) + P g, where U is the particles velocity field, 

= — 1, and P denotes the pressure. The latter is related to the energy density e through an 
equation of state (constitutive relation, conservation law etc0. This is often an ODE between 
e and P, but we will just suppose that P — P{e). 

The compatibility conditions are (ttv is the orthogonal projection to the subspace V) 

{e + P{e))VuU + TTjji^ (gradg P(e)) = 0, Cu{e) + (e + P(e)) dwg{U) = 0, 

and we take these to be the components (ii,L2) of our operator L — L[g,U,€\. The symbol 
^ = CTL : r* (g) (T © M) ^ T e M at ([/, e) has the following wc value on a covector p: 



u 
e 



(e + P(e))(p,C/)./ ^u-{p)P'{e) 



ie + P{e))p 





u 




e 



Since ttj/^ (p) = p — T^jw U we get 



detep = {e + P{e)r-^{p, Uy 



{p,u) 



1 



P'{e) 



This does not vanish for p = C/ and hence does not vanish identically. Thus the operator L is 
determined. Therefore the equations of perfect fluid are involutive. 

2a. Specifications: dust etc. For P = (e > 0) the operator L = [Li^L^) with 

Li = V[7C/-0, L2-diVg([/)+/:t;(loge) = 0. 

has the same form as the one for the null dust. This case is called the dust. Notice though 
that the compatibility condition L = is obtained without re-normalization in this case since 
T = K • f/ ® (contrary to the null dust case). 

Some other cases - incoherent radiation P = e/3 and stiff matter P — e - provide the other 
involutive relativity systems. 

3. Prescribed energy-momentum. Consider now equation ([2|) with fixed non-zero T. We 
change this according to footnotcPl to the equation 



(13) 

(14) 
T at 



G(Ric(5)) = c^*T 

with an unknown local diffeomorphism (p. We add the compatibility (conservation law) 

L[.g,0]=div5(rr) = O. 

It is the 2nd order operator in (1st order in g) and its symbol ~ af^ : S"^ ® T ~ 
(j) € Difrioc(M) equals 

(Tl(p' ® ^){q) = \p\l T{ip, dP{q^)) + (p, q)g T{ip, d^ip^)), p,qeT*,ipe T. 

Clearly if T is degenerate, then <jl{p^ ® ^p) ~ Q for G Ker(T). Let us show that for non- 
degenerate T and \p\g ^ the operator ip = a Lip^ ® ■) :T ^ T is an isomorphism. 

Choose any 7^ 0. If T{(p,d(j>(p^)) ^ 0, then lp{}p){p) = 2\p\lT{(p,d(f>{p^)) ^ and so 
£p((^) 7^ 0. But if T{ip,d(j){pi)) = 0, then ip{ip)iq) = \p\lT{ip,d(l){q^)) ^ for a.e. q e T* since 

: r — > T is an isomorphism. This shows that £p is injective and hence isomorphism. 

Therefore the operator L is determined for non-degenerate T and so the equation £ (fT3|) + (fT4|) 
is involutive. This is a new proof of DeTurck's theorem [DTj . 
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If this is not imposed, the operator L remains underdetermined. 
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6. Conclusion 

In this paper we treated involutivity of several relativity field equations via formal theory of 
PDEs. Some classical cases, such as Dirac-Weyl equations, and many modern specifications are 
not touched, though can be treated via the same technique. Some other cases are more delicate. 

For instance the scalar minimally coupled matter field equations with the action (here ^ is a 
function and dVg — \/ \ det g \ dx the volume form) 

5[.g,0]= / [R-\dcl>\l-m'cj,']dVg 
Jm 

have block-form symbol operator aG®L (with T ~ d4)®d(t)— — |m^(/)^; L[(f)\ — D(j) — iri^ (p) , 

and so arc involutivc by Theorem U 

Non-minimally coupled field equations with the action 

/ [f{cj,)R-\dc^\l~m^cj,']dVg 

have the general symbol operator, and our methods do not apply directly. However the system 
can be re- written with triangular symbol, and then methods of Section [5] show it is involutivj^. 

Let us deduce several corollaries of the involutivity. They are based on the Cartan-Kahler 
theorem claiming that a formally integrable analytic system has local solutions. Since involutivity 
implies formal integrability, we conclude 

Theorem 5. Let j^g he a jet of metric (1 < k < oo), which satisfies [k — 2)-jet of the vacuum 
Einstein equations 0)t=o • Then there exists a local analytic solution g of this equation with the 
prescribed jet j^g at the point £ M . 

In particular, if a Riemann tensor Ricmo at the point is given, which satisfies the obvious 
algebraic compatibility conditions with a metric go G S^TqM through then there exists an 
analytic solution to the vacuum Einstein equations with the given initial data (go,R.iemo). 

This is a variation on Gasqui's theorem |Ga| . For non- vacuum field equation one similarly 
deduce DeTurck theorem jPTj on solvability in analytic category, provided that the tensor T is 
non-degenerate. Turning now to Einstein-Maxwell equation ([9]) or (|10p we arrive at 

Theorem 6. Let j^g^j^F be k-jets of a metric and an analytic 2-form, which are related by 
(jets of) source-free Einstein-Maxwell equation (0),/=o- Then there exists a local analytic solution 
{g,F) of this equation with the prescribed jet {jog, jgE) . 

In particular, if a metric go G S'^TqM , a Riemann tensor Riemo and a 2-form Fo £ A^TqM 
at the point G M are given, which satisfy the algebraic compatibility conditions through the 
first equation of 0), then there exists an analytic solution to the source-free Einstein- Maxwell 
equations with the given initial data (50 , Riemo , -fg ) • 

Finally we can prove local solvability (with prescribed Cauchy data) of other field equations 
considered in Section [5] (pure radiation, perfect fluid, dust etc). 

Theorem 7. Let jog,j'^(j>oF be k-jets of a metric and additional fields, related by (jets of) 
equation il2\) . Suppose that both T and L are analytic. Then there exists a local analytic 
solution {g,(f>) of this equation with the prescribed jet (jo5jJo<^)- 

By similar reasons we have local analytic solutions with any admissible Cauchy data to min- 
imally and non-minimally coupled fields equations etc. 

Appendix A. Some technicalities 

Here we collect some tedious verifications. Let us begin by checking the fact that conservation 
law for ([2]) with fixed T vanishes only for the vacuum case. 

Lemma 1. Let T G S'^T* be a constant tensor. Equation divg(7A) ~ from Section\^ is 
non-trivial only for T = 0. 

Proof. We shall calculate symbol of the operator divg(7A) with respect to g (which varies), and 
regardless T (which stays constant). Metric g enters via Christoffel coefficients, namely in local 
coordinates: 

div,(7X). - (T^l^^ ^ d,^(7A)i + Ti^in)1 - r^(7^)i. 



For almost all functions f{(j>), in particular for the important case /(</i) = 1 — {</i^. 



12 



BORIS KRUGLIKOV 



This yields the following formula for the symbol operator (we g-lift 7a to the operator field) 

This vanishes if and only if the operator 7a is scalar. Since T (as a covariant tensor) is fixed 
((7-independent), this implies 7a = — A g and T = 0. □ 

Next we study exactness of the complexes from Section [2l 

Lemma 2. The operator : S'^^^ ® S'^ S'^ ®T* is epimorphic. 

Proof. We have: 

Since the latter vectors generate S'^ (g) T* the claim follows. □ 
Lemma 3. In the diagrams from Section\^all columns are exact. 
Proof. Formula ([3|) implies that the following sequences are complexes: 

^ 9k+3 3''+^ (g) "^^^ 3''+'^ ® 5*2 S^(g)T* ^ 0. 

Since the previous discussion implies that the only possible cohomology can be supported at the 
term 5'*^+^ (3 S^, the exactness can be verified by dimensional reasons as in |DTi IGa) . 

We will show another approach, but restrict for simplicity of formulas to the case fc = 
(prolongations can be treated similarly). 

The space Ker(crdiv) is generated by vectors p® and \q\^gP®p'^ — 2|p|g q® qp with p ^- q. 
For n> 2 the operator G is invertible: G~^h = h — Trg{h)g. Thus 



Kerf 



(erg) = (1 ® G) ^Ker(CTdiv) = {p®q^ - ^P®g,P'E)q^ - ^p^p"^ + 2q(E)qp\p ±q) 

We need to show that the generating elements belong to Im((T|y|,). This follows from exact 
formulae 



<jlll(^j^ph®pq+ J^2q'^h®ph- j^ph^ - ^^i«iLp2^®p/^+ _|_p3 



•q' 



P®q^ - {^P^P^ + 2(7®gp, 



4l[ 



2 . ^^ „2 

n-2 



{9-p)®q ~ :;^9®{q ®p) + ^^p q®pq + ^^q h®ph 

- i^^Ph^^q^ - '^Wy^^P^ =p®q^ - ^pi 



where p, h are mutually orthogonal elements (with non-zero length; since both kernel and 
image are closed spaces this does not restrict generality), and g®{q^ ® p) = ^£^(9^6^) ® {pci) 
ioT g = J2eief e S^T*. □ 

Finally the following statements was used in Section [S] 

Lemma 4. Consider a vector bundle tt over M with fiber - vector space N. Let F : C°°{t:) 
C°° (M, N) be a determined differential operator of order k. Then the corresponding equation £ — 
Kei{Ap) is involutive and its only non-zero Spencer cohomology groups are H^''^ = H^^^'^ = N . 

Proof. This is a folklore result, and it is difficult to find the reference. One of the proofs is the 
diagram chase by the bi-complex with vertical complex being {g = (Bgt denotes the symbolic 
system corresponding to £) 

0^g^ST*(g)N ^ST*®N^O 

and the horizontal one being the Spencer complex generated by the first column. 

Another possibility is to view this as a direct corollary of Theorem A of IKL4I , which concerns 
reductions of Cohen-Macaulay systems (ideal of the determined equations is clearly Cohen- 
Macaulay). This theorem states that for a subspace V* C T* transversal to the characteristic 
variety Char(£') C FT* (both complexified) the reduced symbolic system g n SV* ® N has the 
same Spencer cohomology as the symbolic system g. It suffices to apply this to any generic 
1-dimensional line V* (for ODEs the claim is obvious) since codimChar(£) = 1. □ 
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